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(57) Abstract: A method and apparatus for use in a hybrid position location system. The method and apparatus combines measure- 
ments from Global Positioning System (GPS) and terrestrial transceiver stations to compute the location of a device. An algebraic 
solution to hybrid position location system equations is output from the method and apparatus. The method and apparatus determines 
the position of a device using a non-iterative method, as against the use of a conventional iterative least mean square method. The 
method of the present invention can be used to solve the location system equations in scenarios where a non-iterative solution is 
desirable. In certain scenarios, the location system equations may have two possible solutions. An iterative method would converge 
on one of the solutions, without any indication of the existence of the other ambiguous solution. Moreover, the iterative method 
may converge on the incorrect of the two ambiguous solutions. Use of the presendy disclosed method and apoarams yjelds^U^e 
ambiguous solutios. The disclosed method mav .be.followed. uP with iterative methods, using the solutions from th^igfe^^etiiod 
Wm^m^^mm^m^m^^^^M^^^od. a different process can then select the correct solution. Thus, the 



algebraic method < 
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METHOD AND APPARATUS FOR DETERMINING AN 
ALGEBRAIC SOLUTION TO GPS TERRESTRIAL HYBRID 
LOCATION SYSTEM EQUATIONS 

5 

FIELD OF THE INVENTION 

The present invention relates generally to locating the position of 
devices, and specifically to a method and apparatus for determining the 
10 position of a device based upon information provided from Global Positioning 
System (GPS) satellites and associated position location systems. 



BACKGROUND OF THE INVENTION 

15 

Recent developments in Global Position System (GPS) and terrestrial 
mobile communications make it desirable to integrate GPS functionality into 
mobile communications devices such as cellular mobile stations. The cellular 
geolocation problem can be solved using either network-based methods or 
2 0 using handset-based methods. 

Terrestrial Location 

Network-based solutions rely on the signal transmitted from the mobile 
station and received at multiple fixed base stations. This can be accomplished 
25 by measuring the Time of Arrival (TO A) of the mobile station signal at the base 
stations. The mobile will lie on a hyperbola defined by the difference in time of 
arrival of the same signal at different base stations. An accurate position 
estimate depends on accurate synchroiuzation and signal structure (bandwidth, 
etc.). 

30 

GPS-based Location 

GPS-based location relies on a constellation of 24 satellites (plus one or 
more in-orbit spares) circling the earth every 12 hours. The satellites are at an 
altitude of 26,000 km. Each satellite transmits two signals: LI (1575.42 MHz) 
35 and L2 (1227.60 MHz). The LI signal is modulated with two Pseudo-random 

1 
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Noise (PN) codes-the protected (P) code and the coarse /acqxiisition (C/A) 
code. The L2 signal ciirries only the P code. Each satellite transmits a unique 
code, allowing the receiver to identify the signals. Civilian navigation receivers 
use only the C/A on the LI frequency. 
5 The idea behind GPS is to use satellites in space as reference points to 

determine location. By accurately measuring the distance from three satellites, 
the receiver "triangulates" its position anywhere on earth. The receiver 
measures distemce by measuring the time required for the signal to travel from 
the satellite to the receiver. However, the problem in measuring the travel time 

10 is to know exactly when the signal left the satellite. To accomplish this, all the 
satellites and the receivers are synchronized in such a way that they generate 
the same code at exactly the same time. Hence, by knowing the time that the 
signal left the satellite, and observing the time it receives the signal based on its 
internal clock, the receiver can determine the travel time of the signal. If the 

15 receiver has an accurate clock synchronized with the GPS satellites, three 
measurements from three satellites are sufficient to determine position in three 
dimensions. Each pseudorange (PR) measurement gives a position on the 
surface of a sphere centered at the corresponding satellite. The GPS satellites 
are placed in a very precise orbit according to the GPS master plan. GPS 

2 0 receivers have a stored "almanac" which indicates where each satellite is in the 

sky at a given time. Groimd stations continuously monitor GPS satellites to 
observe their variation in orbit. Once the satellite position has been measured, 
the information is relayed back to the satellite and the satellite broadcasts these 
minor errors "ephemeris" along with its timing information as part of the 
25 navigation message. 

It is very expensive to have an accurate clock at the GPS receiver. In 
practice, GPS receivers measure time of arrival differences from four satellites 
with respect to its own clock and then solve for both the user's position and the 
clock bias with respect to GPS time. Figure 1 shows four satellites 101, 102, 103, 

3 0 104 and a GPS receiver 105. Measuring time of arrival differences from four 

satellites involves solving a system of four equations with four unknowns given 
the PR measurements and satellite positions (satellite data) as shown in Figure 
1. In other words, due to receiver clock error, the four spheres will not intersect 
at a single point. The receiver then adjusts its clock such that the four spheres 
3 5 intersect at one point. 

Hvbrid Position Location System 

The terrestrial location solution and the GPS solution complement each 

2 
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Other. For example, in rural and suburban areas not too many base stations can 
hear the mobile station, but a GPS receiver can see four or more sateUites. 
Conversely, in dense urban areas and inside buildings, GPS receivers may not 
detect enough satellites. However, the mobile station can see two or more base 
5 stations. The hybrid solution takes advantage of ceUular/PCS information that 
IS abready available to both the mobile station and the network. Combining 
GPS and terrestrial measurements provides substantial improvements in the 
availabihty of the location solution! The hybrid position location system may 
combme Round-trip Delay (RTD) and Pilot phase measurements from the 
1 0 terrestrial network with GPS measurements: 

The hybrid approach merges GPS and network measurements to 
compute the location of the mobile station. The mobile station coUects 
measurements from the GPS constellation and ceUular/PCS network These 
measurements are hised together to produce an estimate of the mobile station 
15 position. 

When enough GPS measurements are available, it is unnecessary to use 
network measurements. However, when there are less than four satellites or in _ 
the case of bad geometry, four or more satellite measurements, 'the 
measurements must be complemented with network measurements Ihe 

2 0 mmimum number of measurements for obtaining a solution will be equal to the 
number of unknowns. Since the system has four unknowns (three coordinates 
and GPS receiver time bias) the minimum number of measurements to obtain a 
solution will be four. For any satellite measurements that are not available 
round trip delay (RTD) measurements may be used to determine the range to a 

25 base station. RTD measurements may also be used to provide time aiding 
mformation. In addition other information, such as PN offset pseudo-ranges (if 
tame bias is the same as for satellites), PN offset differences (if time bias is 
different) and altitude aiding may be used to provide additional information 
and dius mcrease the number of equations that include the unknowns being 

30 sought (i.e.,x,y,z, and time offset). As long as the total number of equations b 
larger than four it wUl be possible to find a solution. 

Round Trip Delav ( RTD) 

The pilot timing on the forward link of each sector in the base station is 
5 ^chronized with GPS system time. TUe mobile station time reference is the 
hme of occurrence, as measured at the mobile station antenna connector of the 
eariiest arriving usable multipath component being used in the demodulation. 
The mobile station time reference is used as the transmit time of the reverse 

3 
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traffic and access channels. 

Figure 2 shows one terrestrial transceiver station 201 and a mobUe 
station 202. As shown in Figure 1, the mobile 202 uses the received time 
reference from the serving base station 201 as its own time reference. 
5 Accounting for its own hardware and software delays, the mobUe station 
transmits its signal such that it is received back at the serving base station 201 
delayed by a total of 2t , assuming that the forward and reverse links have 
essentially equal propagation delays. The total delay is measured at the base 
station by correlating the received signal from the mobile station 202 with the 

10 referenced signal at time T.^,. The measured RTD corresponds to twice the 
distance between the mobUe 202 and the base station 201 (after caUbration of 
base station side hardware delays). 

Note that knowledge of the PN of the serving base station can also be 
used (due to sectorization as a rough angle of arrival (AOA) measurement) to 

1 5 help with resolving ambiguity. 

Pilot Phase Measurements 

The mobile station is continuovisly searching for active and neighboring 
pilots. In the process, it measures the PN offset of each pilot it receives. If the 
20 time reference is the same on both PN offset and satellite measurements then 
the bias on these measurements (as measured at the corresponding antenna 
connector) will be the same. They can then both be regarded as pseudo-ranges. 

If the time references are different then we can simply use PN offset 
differences between each pilot and the reference (earliest arrival) pilot. The 

2 5 pilot PN phase difference is the same as time difference of arrival (TDOA) of the 

two pilots from the two base stations. Figure 4 shows two such base stations 

401 and a mobile station 405. 

Note that on most cellular systems antennas are sectorized and each PN 

is associated with a sector rather than with a base station. Hence, each 

3 0 measurement can provide, in addition to the TDOA information, some level of 

angle of arrival iixformation (AOA) that can be vised to resolve ambiguity. 

Altitude aiding measurement 

It is always possible to determine with which sector the phone is in 
35 commttnicating. This can give an estimate of the phones position to within 
three to five kilometers. Network planning is usually done based on digital 
maps of the coverage area. Based on terrain information and knowledge of the 
sector it is always possible to obtain a good estimate of the user elevation. 

4 
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calculated position, the mobile station 605 reports the pilot phase of the earUest 
arriving path. 

When combinin-; du^ferent types of measurements, iterative solutions 
(such as the weU-known "Nevvton-algorithm" based gradient approach) may be 
5 used to determine the solution (i.e., the position of the device sought). 
However, in certain scenarios in which an iterative solution is used, two 
solutions are possible. Two solutions are possible because of the quadratic 
natiire of the measurements that are used in tiie iterative equation (i.e., the fact 
that at least one of the unknowns for which a solution is required are raised to 

10 ti\e second power). The possible existence of two solutions creates ambiguity in 
the solution. That is, it is not clear which of the two solutions represent tiie 
location sought. This appUes to aU the types of positioning systems (except 
AOA) including the Global Positioning System (GPS). 

The existence of ambiguity dependents on the existence of measurement 

15 redundancy and on ti^e relative locations of the sateUites and terrestrial 
transceiver stations ti\at provide location information. There is always 
ambiguity when there is no redundancy in the measurements. However, 
ambiguity also always exists when there is redundancy, but tiie geometry is 
such that the amount of information provided is insufficient, even in light of 

2 0 additional measurements. However, these are rare occurrences. 

An iterative method will converge to one of \he solutions witiiout any 
indication of tiie existence or position of ti^e other solution. The particular 
solution to which it converges will depend solely on the initial condition used. 

In the case of GPS, because of the distance of the sateUites, the 
25 ambiguous solution is typically very far from the surface of tiie Earth. It is 
therefore impossible that the iterative metiiod would converge to the wrong 
solution if given an iiutial condition close to the surface of the earth. However, 
when combining satellite measurements witia base station measurements it is 
very possible that the two ambiguous solutions will be close to each other. The 

3 0 iterative method would thus converge aibitararily to one of the two solutions 

without a clear determination as to whether the solution to which it converged 
was the correct solution, or whether there are two solutions at all. 

An exhavistive search can be performed to identify both solutions, if two 
solutions exist. However, if only one solution exists, it may be necessary to run 
3 5 the Least Mean Square (LMS) iterative process several times before a 
determination can be made that orJy one solution exists. 

The algebraic method presented by Bancroft ("An Algebraic Solution of 
the GPS equations", published by IEEE on January 8, 1984) and by Schipper 

6 
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SUMMARY OF THE INVENTION 
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Use of the presently disclosed method and apparatus yields both the 
ambiguous solutions. The algebraic method may then be followed up with 
iterative methods, using the solutions from the algebraic method as initial 
estimates of the device location. A different process can then select the correct 
5 solution. Thus, the algebraic method can be used to detect the existence of 
ambiguous solutions, cmd to find both solutior«. 

It should be understood by those skilled in the art that the method and 
apparatus disclosed is described in the context of a hybrid GPS and ceUular 
location system. However, the disclosed method and apparatus is equally 

10 applicable to any location system that combines satellite and terrestrial 
measurements, such as integrated GPS and Long Range Navigation (LORAN) 
or other such terrestrial systems. 

The present invention will be more fully imderstood from the following 
detailed description of the preferred embodiments thereof, taken together with 

15 the following drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

2 0 Figure 1 shows four satellites and a GPS receiver; 

Figure 2 shows on terrestrial transceiver station and a mobile station; 

Figure 3 shows three satellites a terrestrial transceiver station, and a 
mobile station; 

Figure 4 shows two such base station and a mobile station; 
25 Figure 5 shows two satellites 501, 502, a base station 504, and a mobile 

station 505; 

Figure 6 illustrates a satellite 601, two terrestrial transceiver stations 604, 
and a mobile station; and 

Figure 7 shows the structvire of one example of a device used to 

3 0 implement the disclosed method and apparatus. 



DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

35 Overview 

The disclosed method and apparatus is a system that uses both terrestrial 
transceiver stations and satellites (i.e., a hybrid position location system) to 

8 
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system m which there are either not enough satellite measurements to 
5 de tr^V " " ^ --rate position^ ^ 

5 ^"-"usmgacombinaHonofsatemtesandterr^trialtranLiverstTor 
such as base staHons of a ceUular communication system ' 

In accordance with the disclosed method and apparatus, an "aleebraic" 
method . used to determine whether two solutions ex^^and ttteLue^'ll 

preferable for obtammg both ambiguous solutions. The disclosed method and 

^!ZT< - approximate, Clutn tTa 

system of navgahon equations. system of navigation equations incudes 
one equahon for each of the foUowing: (1, the altitude of thT^ 
d«ermi„ed by altitude aiding information. 2 s.^m^Zl:^^t,^';^ 
' -dmg ^formation (i.e.. an estimate of the receiver docK bias); ~te™ 

sv^temTr ^'^'^^^ - ^ used to™ 

' " -"^'^ a «.I"«on is 

^'^^ 'PP"'^'«°'' proposed here reUes on linearizing satellite and 

.ranL^~" T""*"^ °' --ori^d terrestrial 

E911 pW rZ " unicating (equivalent to informaHon for 

other mt^ or^":^ ^ """^ ^^^^'^ 

fixes, info^aHon JZTd'tfr " ^"^'^ 

another embod";;' T 

is either the cr^offte^ ^«">-« 

Accordingly T^^d ^ T"^ T" " """"S 

mgiy, .t should be understood that the estimate can be m»H» ■ 

mformatron regarding the location of the serving secteTa^^ ^ 
StaLTe^our -"'^ P-n rrbt^r, 

^icaurb:':;rte~^^tr°"^'"--^«^--- 

linearizing the sateUite and altitude aiZTn, Wro=amation made by 

^Sebraic location deter^tion :::^tr^Z:.7Z'l^Z 



wo 01/48506 



PCTAJSOO/33375 



10 



that are raised to a pov/er of two (i.e., unknowns of the second order) can be 
grouped together to form a single variable. The variable must be defined the 
same way in each of the navigation equations. This is not possible in the case in 
which the four types of navigation equations noted above are presented due to 
the differences in the form of each of these four equations. Linearizing the 
sateUite and altitude aiding measurements reduces the number of second order 
unknowns and thus allows the second order unknowns to be grouped together 
and defined as a quadratic variable having a consistent definition throughout 
each of the navigation equations. 

If one of the solutions is more than 15km from the reference point, then 
that solution will be inaccurate. However, such an inaccurate solution will not 
be the desired solution, since we have predetermined that the solution should 
be within 15km of the reference point. In cases in which the reference point 
cannot be predetermined to be within 15km of the user (i.e., in systems in which 
1 5 cells have a radius of greater than 15 km, such as in AustraKa), the accuracy of 
the approximation can be improved if the altitude information is not 
approximated by a plane wave approximation. 

Accordingly, if the altitude information is linearized, then the ambiguity 
can be resolved as long as only one of the ambiguous solutions is within 10- 
2 0 15km of the center of the reference point. If both solutions are within 15km of 
the reference point, then the approximations are vaUd for both solutions. 
Therefore, the estimates for both solutions are accurate and one solution cannot 
^be sf l(^ed over the other. Therefore, other criteria must be used to distinguish 
die desired solution from die erroneous solution. 
25 Once an approximate solution is determined, the approximate solution 

can be used as the imtial condition for determining a more accurate iterative 
solution. Using the solutions from the approximate solution as irutial estimates 
of the mobile station location provides a rapid convergence to a solution that 
lacks the error introduced by the approximation. 
30 Some of the criteria that can be used to identify the correct solution 

include, but are not limited to: (1) sector angle opening (i.e., the angular size of 
the sector) and orientation, (2) distance to serving base station relative to 
expected ceU size, (3) relative LMS cost of the two solutions in the case where 
there is redundancy, (4) received signal power and (5) Coverage maps available 
3 5 for network planning. Coverage maps would constitute the optimal criterion. 

Although the description of the method in this document uses a hybrid 
GPS and cellular location system as an example, it can easily be applied to any 
location system that combines satellite and terrestrial measurements, such as 

10 
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integrated GPS and LORAN. 
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navigation equations assume that an estimate of the receiver position is 
available that is accurate to within 10-15kms. In general, the sector with the 
earUest time of arrival at the phone will be referred to as the serving sector. The 
referertce point wiU be the center of the coverage area of the serving sector. 
5 Note that if the size of the sector is larger than 10-15 km, then it may be 
necessary to run an iteration of the disclosed method wherein the reference 
point is updated according to the result. In general however this wUl not be 
necessary. 

10 Altitude Aiding 

An estimate of the altitiade of the mobUe station may be available from 
terrain information, previous location solutions, or other sources or 
measurements. If the mobile location x„=[x„ y„ zj is defined in Earth 
Centered Earth Fixed (ECEF) coordinates, the estimate of tiie altitiide is an 

15 estimate of |3c„| . In order to include altitiide aiding in ihe algebraic solution, we 
must express tiie altitiide aiding equation as a linear equation so that it will not 
resti-ict the choice of the algebraic method quadratic term. This can be achieved- 
by rotating the coordinate frame such that the altitiide estimate becomes a 
linear combination of the unknowns in the system of equations (witiiin a certain 

2 0 radius of the rotation reference point). 

We rotate the ECEF coordinate frame such that the Z-axis passes through 
a point selecf d as an initial estimate of the location of the mobUe station. In the 
case of a hybrid location system using terrestrial tiarisceiver statior«, such as 
cellular base stations, and^lSfTikteUites, this initial estimate can be a point in 
25 the coverage area of a selected base station. If the base station measurements 
are pseudo-ranges or ranges then the center of the serving sector can be used as 
initial estimate. If the base station measurements are range-differences, then the 
serving base station (range-difference reference) will have to be used as the 
initial estimate. This is due to the constraints imposed by the method in the 

3 0 case of range-difference measurements. 

If the irutial estimate of the mobile station location is dose to the true 
location of the mobile station, then the estimate of the mobile station altitude is 
an estimate of the Z coordinate of the mobile station in the new rotated 
coordinate frame. It is obvious that linearization can alternatively be 
35 accomplished by transforming the altitude estimate to an estimate of the X 
coordinate or the Y coordinate of the mobile station (instead of the Z 
coordinate, as described above). The rotation matrix T is computed as follows: 
If ^0 =[;Co yo Zq] represents the ECEF coordinates of the irutial estimate 
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SDheriral rn«r^;«,.. ... « " « • <P and r are the coordmates in the 



spherical coordinate frame, then 

(3) 



7- = 



cos(^)cos(^) cos(^)siii((!)) -sin(^) 
-sin($») cos((*) 0 

Lsin(^)cos(«>) sin(6')sin(<!>) cos(^) 



(4) 



? represent. Ae coordinate, of satellite ; in ECEF coordinates, and I 
represents the coordinates of satellite i in the r<,b,t^ 

coordinates/^n^rt^tr ^rvirdTrr " T ""'"^ " 
station in the rotated coordinate sCZ 

be taken into account sZTbvTdd °' ^ """""'"^ ^ 

system of equado^ Zl^' of t^' ' ^ 
equations in\r:ns of T Ef^tiot ^ " '° ''^'^^ '^^^ 

selected to make the equality true. '- 



^ a a 



/. is equal to zero, and c. is equal to i . unknowns x, y, z, and 6, 
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Au=l X + c <=>[0 0 1 0 

A a a 



= 0»X + z - (6b) 



The form of Eqiiation 6b makes it easier to combine the altitude 
information with the other information that is known, such as the satelUte 
measturements and the base station measurements, as will be seen below. 

5 

Satellite Measurements 

Let [x„ y„ z„] be the location of the mobile station whose location is 
sought and let [x„ zj be the location of a satellite S, . Let b be the receiver 
clock bias. Accordingly, the pseudo-range measurement to each satellite, 
10 p , where i = 1... ., n , can be expressed as: 



P„- =a/U, J +fc -yJ '■=''•••'' 

Because the satellites are far away from the Earth, it is reasonable to use a 
15 plane-wave approximation. The plane-wave approximation assumes that, 
instead of a sphere at a distance from the satelUte, the sateUite measurement 
surface is a plane at a distance from the satellite. 

A vector is defined as v,, = ,l^~l" , /the line of sight vector going from 

\Xr-X,i\ 

the satellite to the reference. The satellite measurement equation can be written 
2 0 as: 



<x-x ,v >+b = p (8a) 

JI si si 

It will be understood that Equation 8b follows from Equation 8a. 
<jc,v >+fr = p +<x ,v > (8b) 

51 si si si 



30 
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A system of equations can be written in fKo f^u • r 
the relationship of each set of sateZ^ "^"^ 
sought: measurements to the location that is 



(9a) 



10 



v^.Cl) v;_(2) v^(3) 1 
v7^(2) 57^^ (3) 1 

\„(2) i7 (3) 1 



P+<x ,v > 

*2 j2' ,2 



P +<a: ,v > 

" mm 



(9b) 



Time AiHing 

es«n.at::.x~"::eror - "> 

measuring the amount of dme reaJLed f T made by 

station ,o reach a mobUe sta^I b ' ' transmitted f„,m a base 

received by ,he bCe station ^^""""'^ *e mobile station, and be 

transmined and ™ "Lt a"' retransmission (i.e.. *e 

assumpHon is ma^Zt tofT ^'™^™°>«> mobile station. If an 

the amount of^e r^'fJ^'T " ""^ 

mobUe station can bTd^rC hLr.^;^ 

^rztf:^«:e^;r r° " 

base station and rmX l '° *^ "^'""^ 

base staHon can be^e:™ ;Ttdr°" ^ "> 

is used as the time referrcHo m t ^'"^ 

tune reference to measure the GPS pseudoranges. Accordingly 



RTD 
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be available from other sources or measurements, such as previous navigation 
solutior\s. The estimate of the clock bias can be taken into account simply by 
expressing a new linear equation, fe = i?, in the form used for the system of 



equations. Here, Aj-u ^l^A-^c^ <;=>[0 0 0 1 



= OxA + t (11) 



Terrestrial Measurements and System Resolution 

The terrestrial measurements can be treated in three v^ays: 

1. Terrestrial pseudoranges 

2. Terrestrial ranges 

10 3. Terrestrial time difference of arrivals. 

Terrestrial measurements as pseudo-ranges 

Pilot phase measurements made by the mobile station can be treated as 
pseudoranges. In a system using both GPS and LORAN, LORAN 
15 measurements may be treated as pseudoranges. If the terrestrial measurements 
are treated as pseudor£uiges, they can be expressed as: 

2 0 v^here is the bias on each measurement. 

For each one of the measurements we perform the following 
manipulation. First subtract b, from both sides. Then square both sides of the 
equation, resulting in: 

25 



(p -b}^\x'X r (13a) 
^bi \ bi\ 

Next, each side is expanded, resialting in: 

< - 2/^.,^, ^ K = 1^' - 2 < X, J > H^^J' (13b) 

Next, all of the second order terms are collected on the right side of the 
4 0 equation: 

16 
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2<x,x 



i,>-^.*=H'-*=%J-/t. (130) 



A quadratic variable, A=|l|'-i,=, is defined. Equati, 
written as: 



on (13c) can be 



2<x, 



25 



We can perform the operations in Equations (13) and 
measuremei 
ons for a plu 
in the desired form as: 



(14) 

(14) on the pseudo- 



2y,, 2z,, -2p,2 
.2^,„ 2y,^ 2z^ -2p,„_ 



1 



(15) 



defined^^H^rHoUrS)!,!^^^^^ " - 

at this stage Noru^a 1 j V""''"'™'"' """^'^ 
approximaHonfor ftralTrf^ " """^"^ '° Plane-wave 















A, 


- h 

u~ * 
0 








A. 


[o_ 


1 





(16) 



Let B be the generalized inverse of A rnnt^ rt,=.f *k 
this case is not the same «c f>,o ^ covariance matrix in 

the same as the covanance matrbc of the measurements), then: 



30 



'*-Bx(lX + c)=BxU + Bxc = pX^q 



(17) 
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10 



20 



25 



We define a par of vectors d and e that represent the x, y, and z, 
components of the vectc rs p and q as follows: 

e=(9(l) ^(2) ^(3)f 

and a pair of scalars/ and g that represent the offset h of the vectors p and 
q as follows: 



30 



^/ = P(4) 
1^ = 9(4) 



(18b) 



Therefore, we can see that: 

x = dX + e (18c) 

15 p|'=p;i + ;|' (18d) 

Therefore, substituting Equation 28b into the definition of X results in: 
A = |3cl'-fe;=pA + 6l'-(/A+g)^ = 



(19a) 

l^'A^+2<J,e>A + |?|'-(/U^+2/^A + gO 



We then collect all of the terms associated with X* together, all of the 
terms associated with X, and all of the terms that are imassociated with X on the 
left side of the equation. 

\^d^ - + (2 < J,e > -2/5 - + 1^1' - = 0 (19b) 

Equation (19b) is a second order equation in X cind has the following 
solutions: 
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^_ -{2<d,e> -2fg - 1)+ j(2 < > -2/8 -^f -^d f - f^^' - g^) 

^ _ -(2<J,e> -2fg - 1)-^(2 < J,e > -2/g - if -^^f^^^ f^'^) ^""^ 

We can find the solutions that correspond to these two values of J, by 
^ ^-bstitutmgthemintothedefinitionofthesystemvariablesasafunctionofX: 

u =pA +q (21) 

back into If ^■»"«°"- ''"b^tttute ftese hvo solutions 

I^hlr^' '"^''^ *^ ^y^t^n has two 

ambiguous solutions. 

Terrestrial measnrPT nents as ran |arpc 

*e distance between the mobile slaHon and the base staHon. Ihe RTD 
— ement made at the refere„« base station can be treated as a range 
measurement T„e RTD measurement made at the reference base slahon can be 
oombmed with mobUe station's measurement of the «me difference of a3 of 

a.titude'rilXiraaWe'^' " — - 

expresl*V:~ "^^^^^ "^^'^ ^ - 



25 



30 



manipulation 



(22) 

For^^each one of the measurements we can perfonn the following 
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bi \ hi 



(23a) 



10 



15 



20 



The term on the right side of the equation is then expanded: 
r^=l3c|'-2<*,jc >-»-|3c T (23b) 

bi I ' bi \ bi\ 

The second order terms are then collected and isolated on the right side 
of the equation: 



2<jc,jc >=l5cl 
bi ' ' 



X -r 

bi I bi 



(23c) 



A quadratic variable is defined as A = |jc|\ Equation (23c) can be 
expressed as: 



(24) 



We can perform the operations in Equatior\s (23) and (24) for all the 
range measurements. The system of equations can then be written as: 





2>M 


2Zm 


0' 


X 




"f 










2>'fc2 


2Zi2 


0 


y 




1 


A + 


\Xb2 












z 




1 










2^*, 


22b„ 


0 


b 




1 









(25) 
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The altitude aiding, satellite and time bias measurement equations as 
defined in Equatioiis (6), (9), and (11), respectively, can be added to the system 
at this stage. The sets of equations can be concatenated so as to obtain a single 
set of equations: 





As 




0 






Au=lX-\rC = 


A, 


u = 


h 








At 




0 








A. 




0 







(26) 



Let B be the generalized inverse of A (note that the covariance matrix in 
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this case is not the same as the covariance matrix of the measurements), then: 
"=Bx(U + c)=BxU^Bxc = pA + q (27) 

We define a pair of vectors d and . that represent the y, and z 
components of the vectors p and g as foUows: 

R = [p(l) p(2) p(3)Y 

[e=[m g(2) ^(3)f ' (28a) 

of ^ ^^"^ °^ '"^'"^ ^ ^^P^"^^"* o«^et & components 

of the vectors p and as follows: omponents 

f/ = P(4) 

U=^(4) (28b) 
15 Therefore, v^re can see that: - 



10 



20 



25 



|2 -12 



(28c) 

H'=R^-»-? (28d) 
Therefore: 



zero: 



^-m''J^A*ef^P\'A'^2<S.e>i^^' (29a) 
By subtracting X trom both sides of Equation 29a, the^quaUty is set ,o 

Pl'-l'+(2<d,r>-i>l+|y|'=0 (29^ 



° soluHol'"^*" ^ in . and has Mowing 
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-(2 < d,e > -l)+^(2 <d,e> -\f -4|J|^|«|^ 



-{2<d,e>-l)-^(2<d,e>-lf -^df\e\^ 



(30) 



We can find the solutions that correspond to these two values of X by 
substituting them into the definition of the system variables as a function of X: 



u = p X -h q 



(31) 



10 



To distinguish the correct solution, we substitute these two solutions 
back into the system of equations to find the solution that yields very small 
residuals. If both solutioris yield small residuals, the system has two 
ambiguous solutions. 



Base station measurements as range differences 

Hie mobile station measures the time difference of arrival of the pilot 
15 signal from different base stations. These measurements can be treated as range 
differences. In a system using both GPS and LORAN, LORAN measurements 
may be treated as range differences. We assume, without loss of generality, that 
one of the base stations (say b^) is the reference for all the range difference 

measurements and that this base station is the origin of the coordinate frame. 
2 0 Therefore, the range difference measurements can be expressed as: 



(32) 



25 



For each one of the measurements we perform the following 
manipulation: 



(33a) 



Expand both sides: 
Sl+2S\x\-^M' -2<x,x >+\x 
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ion: 



Then coUect the second order tenns on the right side of the equatic 
2 < J, X >= -2 J Ijcl + |x r (33c) 

bt* * \ bi\ bi ' 

A quadratic variable is defined as A = |x| . Equation (33c) can be written 



10 as: 



2<x,x,,>=-2S,,A + \x,,\'-S^, 



(34) 



We can perform the manipulations in Equations (33) and (34) on all the 
range difference measurements. The system of equations can then be written 



as: 



20 



^^bi 2yt,i 2zj, O" 
2*fc2 2>„ 2z,2 0 

2x,„ 2y,„ 2z,„ 0 



-2S. 



b2 



-2S, 



in J 



— |2 
^b2 



-SI, 



(35) 



The altitude aiding, satellite and time bias measurement equations as 
defined in Equations (6), (9), and (11) can be added to the system at this stage. 
The sets of equations can be concatenated so as to obtain a single set of 
equations: 



25 
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Am = / A + c = 







"o" 






As 


u = 


lb 
0 






A. 




0 







(36) 



Let B be the generalized inverse of A (Note that the covariance matrix in 
this case IS not the same as the covariance matrbc of the measurements), then: 



(37) 



23 

BNSDOCID: <WO 0148506A2_I_> 



wo 01/48506 PCT/USOO/33375 

At this point we define two new vectors d and e and two new scalars / 

and g: 

;j'\pm mr ^jlz-fw (38) 

P(0 P(2) p(3) y gj^^ 2 components of the vector p. 

q (1) q (2) ^ (3) £ire the x, y and z components of the vector q, 

p{A) is the h component of the vector p. 

q{A) is the h comj>onent of the vector q. 

This allows the x, y, and z components to be treated separately from the b 
component. 

If we substitute Equation 38 into Equation 37, we get: 

M = + ^ = J>1 + e + /A -f g 
It should be seen that: 

x = dX-¥€ 

Accordingly, the system imknovms x, y, and z can be expressed as a 
function of X as follows: 

A = jjcj » A' = |3c|' = pA + = A' + 2 < > A + |el' (39a) 
Equation 39a is set equal to zero by subtracting X from both sides: 
\^ - 1 j^^ + 2 < d , ^ > A + = 0 (39b) 
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Equation (39b) is a second order equation in X, and has the following 
solutions: 



^ _ d,e>-p<d,e >f - 4^\df -l^f 



(40) 

2<. 



We can find the solutions that correspond to these two values of k by 
substituting them into the definition of the system variables as a function of X: 

u '^pA +q (41) 

2 2 ^ 



Ambi^ty Resolutinn 

Algebraic resolutions of quadratic systems always yield two solutions, 
even in case of redundancy. In accordance with one embodiment of the 
disclosed method and apparatus, to distinguish the correct solution, we 
substitute these two solutions back into the system of equations to find the 
solution that yields small residuals. If both solutions yield small residuals, the 
system has two ambiguous solutions. The correct solution wiU be the one that 
is consistent with the sector information associated with the base station 
measurements. Alternatively, it will be understood by those skilled in the art 
20 that any of the method used to determine the initial estimate of the location 
sought may also be used to assist in resolving the ambiguity (i.e., select one of 
the two solutions). For example, the sector that is in communication with the 
device whose location is being sought may eliminate one of the solutions, 
alternatively, the location of the serving base station, the altitude of the device 
as determined by an altitude sensor within the device, or any other information 
that might be used to limit the possibility that one of the solutions is more Ukely 
to be correct. As noted above, some of the criteria that can be used to resolve 
the ambiguity include, but are not limited to: (1) sector angle opening (i.e., the 
angular size of the sector) and orientation, (2) distance to serving base station 
relative to expected cell size, (3) relative LMS cost of the two solutions in the 
case where there is redundancy, (4) received signal power and (5) Coverage 
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maps available for netiA ork planning. 

Figure 7 shows the structure of one device 700 used to implement the 
disclosed method and apparatus. As shown in Figure 7, the device 700 includes 
an antenna 702, a transceiver 704, and a processor 706. The antenna receives 
5 signals from each of the signal sources, such as satellites and terrestrial 
transceiver stations. The signals are coupled from the anterma 702 to the 
transceiver 704. The signals are then processed by the transceiver 704 in a 
manner well-known to those skilled in the art. The transceiver may be an 
analog communications transceiver, digital communications transceiver, GPS 
10 position location transceiver, Loran transceiver, or any combination of these or 
other types of transceivers. The processed signals are then coupled to the 
processor 706. The processor 706 may be any type of computational device that 
is capable of performing the functions described above, including a general 
^ purpose microprocessor including memory, a special purpose microprocessor 
15 including memory, an application specific integrated circuit (ASIC) (or portion 
of an ASIC), dedicated circuitry comprising discrete components, a state 
machine, or any general purpose computer, including mini-computer, desktop 
computer, laptop computer, or mainframe computer. The processor 706 
outputs the location of the device 700. It should be xmderstood that the 

2 0 processing functior\s performed by the processor 700 may be distributed among 

several components that may or may not reside in the same physical location. 
For example, it is common for information to be collected by a device and 
transmitted to an external device, such as position determination equipment 
(PDE) which performs some of the required calculations and marupulations. 
25 It should be noted that the preferred embodiments described above are 

cited by way of example, and the full scope of the invention is limited only by 
the claims. For example, while the application notes the use of communication 
base stations in several examples above, the terrestrial transceiver stations may 
be any station capable of providing signals that would accommodate the 

3 0 current method and apparatus for determirung position location. Likewise, the 

satellites referred to in many of the above examples are GPS satellites. 
Nonetheless, it will be imderstood that the satellites may be any system which 
provides additional signals that provide position location that can be used as 
described above for position location determinations, 

35 

WE CLAIM: 
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CLAIMS 

1. A method for determining the location of a device including the steps of: 

2 a) receiving range information, pseudo-range information, and difference 

of arrival information related to a particular location sought to be 
4 determined; 

b) using a plane-wave approximation to eliminate the xmknown second 
6 order terms associated with the pseudo-range information; 

c) substituting a quadratic variable for the xmknown second order terms in 
8 the range information; 

d) constructing a coordinate frame with one of the transmission points as 
10 the origin of the coordinate frame; 

e) expressing the range difference of arrival information in terms of the 
12 newly constructed coordinate frame; 

f) substituting a quadratic variable for the coordinates of the unknown 
14 location, thus placing the equation for the difference" of arrival 

information in the same form as the form used to express the pseudo- 
16 range information and the range information; 

g) concatenating the equations for the range, pseudo-range and difference 
18 of arrival information into a single set of equations; 

h) expressing the unknowns as a function of the quadratic variable; and 

20 i) solving for the quadratic variable, and thus determining two solutions 

for the location sought. 



■ <WO_0148506A2J_> 



27 



WO Ol/48506 ^ ^ ^ PCT/USOO/33375 




HGUREl 



BNSDOCID- <WO 0148506A2_I_> 



wo 01/48506 PCTAJSOO/33375 

2 / 7 





BNSDCXJID: <WO 0148506A2_L> 



PCT/USOO/33375 

WO 01/48506 3 /7 




BNSDOCtD- <WO 0148506A2_I_> 



wo 01/48506 



4 / 7 



PCT/USOO/33375 




Base Station 1 



FIGURE 4 



: <WO 0148506A2_I_> 



PCT/USOO/33375 

WO 01/48506 . . 




505 



Mobile Station 

FIGURE 5 



BNSDOCID <WO 0148506A2_I_> 



wo 01/48506 



6 / 7 



PCT/USOO/33375 




HGURE6 



BNSDOCID: <WO 0148ErdA2_l_> 



wo 01/48506 



7 / 7 



PCTAJSOO/33375 




BNSDOCID <WO 0148506A2_L> 



(USPTO) 



,„„«ERNAT,0NA.^UCAT.0NPUBUSHO.«m.Ej.THEPATENTC0OPEKA„0«TREATV(PCT» 



(19) World Intellectual Property Organization 

International Bureau 

(43) International Publication Date 
5 July 2001 (05.07.2001) 




PCX 



(10) International Publication Number 

WO 01/48506 A3 



(SI) international Patent Classifcation-. GDIS 5/14, 5/12 (74) 
(21) International AppUcatlon Number: PCT/USOO/33375 



(22) InternationalFilingDate: 7December2000(07.l2.2000) 

English 

(25) Filing Language: 

(26) Publication Language: 



(81) 



English 



10 Decemberl999 (10.12.1999) 



US 



(30) Priority Data: 

09/460,180 

(71) Applicant: QUALCOMM INCORPORATED [US/USlV 
5775 Morehouse Drive. San Diego, CA 92121-1714 (US). 

Inventors: FERNANDEZ-CORBATON, Ivan, J.; 919 

oTamtd Street, San Diego. CA 92109 (US). VAYANOS^ 
Alkinoos, Hector; 1134 South Fels^ Street #2 San 
Diego. CA 92109 (US). AGASHE, 10173 
Cat^iio Ruiz #4, San Diego. CA 92 126 (US). SOLIMAN 
Samir, S.; 1 1412 Cypress Canyon Park Dnve. San Diego, 
CA 92131 (US). 



(72) 



Agents: WADSWORTH. Philip, R. et al.; Q-^lco^ In- 
corporated, 5775 Morehouse Drive. San Diego. CA 92121- 
1714 (US). 

Designated States (national): AE, AG. AL AM.^. AU. 
AZ. BA, BB. BG, BR. BY, BZ, CA. CH. CN. CR. CU. CZ. 
DE DK, DM, DZ. EE. ES. Fl. GB. GD. GE GH GM, HR, 
HI I tn IL IN IS, JP, KE. KG. KP. KR. KZ. LC, LK, LR, 
S^'lT LU IV MA, MD, MG, MK. MN, MW, MX. MZ, 
n6 i ft. RO. RU. SD. SE, SG, SI. SK, SL. TJ. TM. 
TrJiT, TZ, UA, UG. UZ, VN. YU, ZA, ZW. 

at\ Designated States (regional): ARIPO patent (GH. GM, 
^ ^ S MW. MZ. SD SL, SZ. TZ, UG, ZW). Eurasian 
Snt(^.AZ.BY,KG,KZ.MDRUJJ/™XEuro^^ 
Satent (AT, BE. CH. CY. DE. DK, ES H, FR GB. GR. ffi. 
m LU. MC. NL, FT. SE, TR), OAPI P^«« (B^^i:^:^- 
CG, a, CM. GA. GN, GW, ML, MR, NE. SN, TD, TG). 

Published: 

with international search report 

(88) Date of publication of the international ^^^-^^^^^ 



[Continued on next page] 



;; .,T,a.- METHOD «^.PP^^SH>..D.n.»M»a,« AN « ^»A.CSOLm,OH^C^^RRBStK^^ HV- 
BRID LOCATION SYSTEM EQUATIONS 



5 
o 

ID 



o 




(S7)Abstract:An«thoaandapparatusforuscinahybn.„^^^^^^^^^ 

2t; from Global Positioning System (GPS) and ^--t^^^^l^ST^^^^ method and apparatus d^es 

solution to hybrid positionlocation system equauoi^isoutprnto 

Ae position of a device using a non-iterat.ve method, ^J^"^ ' ^^ons in scenarios where a non-iterauve soludon is 

metLd of the present invention can be used to solve ^^^f ^^^^^^^^^ solutions. An iteraUve method would convey 
desirable. In certain scenarios, the location ^^>^|f "L^^^^.^^C^^^^^^ solution. Moreover, the itera-.^ ™e^h^^ 

on one of the solutions, without any mdication °f "^^^^ ^le presenUy disclosed method and apparaws yields both the 
may converge on the inconect of the two ambiguous soluaons^^ 

amLguouslolutios. m ^^-^^-^ -^^^JJ,*;,^^^^^ select the correct soluUon. Thus, the 

ssSctrrait^'s™*^^^^^^ 




For two-letter codes and other abbr^iations. refer to the "Guid- 

Zl f"'T ^1"' ""^^''''^""'ons-appearingat the begir,- 
ning of each regular issue of the PCT Gazette. 



BNSDCCID: <WO_014850eA3_L> 



INTERNATIONAL SEARCH REPORT 



int Hional Application No 

PCl/US 00/33375 



.n.P.en.aos.lca..onOPC a, c.ass,..3Uon a.. . PC. 

IPC 7 GOIS 



Docu menial ion 



ITarched other than minimum documemauun lo 



The exient that such documenis are 



included m ihe neids searched 



EPO-Internal . WPI Data, INSPEC 



"where praciicai. search lerms useo) 



r nOCUMENTSCOr^'n^HEDTOBEREl-EVANT 



Category " 



UO 99 54752 A (SNAPTRACK INC) 
28 October 1999 (1999-10-28) - 

page 19, paragraph 3 -page par.»y 
BANCROFT S: "AN ALaiBRAIC SOLUTION OF THE 

us, 

ofnuaffllsB 0986:01), pages 56-69. 

XP000972710 
ISSN: 0018-9251 
cited in the application 
the whole document 



I Qj F„rtherdocun,en.sare.B.edinlhe con..nua.K,n ot box C 

P Special categones ol cited documenis 

A- documenldetmlng the general slale 01 the anwh^h.shol 

^SS ed 10 he Ol particular relevance 
•E- earteTL.a«n.bu.puh.lshedonora.,erlhe ,n,ema.iona. 

•O- ru:^nUe.ernnro-or^<«^-» "- 

•p. ^^^hl^edPyr^^^^^^^ — 
later than itie pno niy dale claimeo _ 

[-5a,e 01 me ac.i.a. coniple.ion ot the ,n.erna..«..-i sea.i.1. 

12 November 2001 

--'"-•™Trot^Pa.e":.'S:.e.P.B. 5Bte Palenllaan . 

reW.?t-oT3.«O.Tx3t 
Fax; (+31-70) 340-3016 



-/— 

ID 



Patent lamily members are Isted in annex. 



later documen, Pu-^'^ ^'Ji'S v!S^^ 
a^u^rrst^d m°LS^^^ ""-^^-^ - 

. <^":::^"c.pad^--^3"-^^^^^^ 

: d^cu-nlotparti^larreleva^^^^^^^^^^^ -"rnthe 
cannol be considered to involve an " ^jq^u- 
rn.T?.:^Srinrn'X^°=.oa person sailed 

I- documeni member ot ihe same patent famity 
Date ot ma.i.ngot the mternatwnat search repon 

20/11/2001 

Auinonzed otttcer 



Roost, J 




:ofTT, PCT/lSA/210 (second shoet) (Juiy 19921 



page 



BNSDOCID- <WO 0148506A3J_> 



INTERNATIONAL SEARCH REPORT 




1 

Fom> PCT/ISA«10 (conimual.on m ueona siwelKJuly laszi 
BNSDOCID: <WO 0148506A3.L. 



page 2 of 2 




Form PCT/rSA/210 (patent family anne«HJuly 1992) 



